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Введение. Предложена концепция оценки динамических параметров системы «основание — ослабленный слой — 
блок» с учетом физической нелинейности материала при кинематическом способе возбуждения колебаний. В 
соответствии с данным подходом учет физической нелинейности материала основания и блока осуществляется с 
помощью модели Друкера-Прагера. Ослабленный слой моделируется ЗО пружинными конечными элементами. На 
примере динамического расчета системы «основание — ослабленный слой — склон» осуществлена процедура 
верификации предлагаемой методики. 

Материалы и методы. Вычислительные эксперименты выполнены с помощью программного комплекса АМЗУЗ 
Месватса| в сочетании с нелинейным решателем, базирующемся на процедуре Ньютона-Рафсона. Для 
дискретизации расчетных областей применены объемные конечные элементы $045. Для моделирования 
смещения блока относительно неподвижного основания использованы комбинированные упруго-вязкие элементы 
СОМВИ\ 14. 

Результаты исследования. Разработана инженерная методика динамического анализа напряженно- 
деформированного состояния пространственной системы «основание — ослабленный слой — блок» при 
кинематическом способе возбуждения колебаний. На числовых примерах исследована точность и сходимость 
предлагаемой методики. 

Обсуждение и заключения. На основании выполненного математического моделирования показано, что 
разработанная методика позволяет оценить риски возникновения реальных оползневых процессов, обусловленных 
внешними нестационарными воздействиями. 

Ключевые слова: метод конечных элементов, модель Друкера-Прагера, метод Ньютона-Рафсона, динамическое 
моделирование, склон, оползневой процесс. 
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Введение. В 70-х годах прошлого века началось широкое использование метода конечных 
элементов (МКЭ) при анализе статической устойчивости склонов и откосов [1-3]. Как показала практика 
геотехнических расчетов МКЭ в отличие от упрощенных методик позволил учесть такие важные факторы как 
реальную геометрию и послойную структуру исследуемых объектов, а также наличие противооползневых 
сооружений и физическую нелинейность материала с заранее неизвестной областью пластичности [4]. В настоящее 
время многие исследователи стали применять МКЭ для анализа динамической устойчивости склонов реального 
заложения, а также прогнозирования рисков обвалов горных пород в подземных выработках и со склонов, 
расположенных вдоль автомобильных и железнодорожных магистралей. Также одним из актуальных направлений 
горной динамики является конечно-элементное моделирование последствий землетрясений с учетом 
кинематической нестабильности конгломераций [5]. В настоящее время существуют следующие методы 
динамических расчетов в области геотехники: 

1. ВМ (зепо@ тедисйоп тео) — метод снижения прочности [6, 7]. Разработан для расчета запаса 
прочности горного массива в физически нелинейной постановке. В ЗКМ фактические прочностные параметры 


грунта итерационно уменьшаются путем деления на некоторый коэффициент больше 1: 








С ап 
Е = ; ф’=агсапт р | 
КИ 5АМ 
где с, ф — фактические значения сцепления и угла внутреннего трения грунта соответственно; с’, ф’ — 


сцепление и угол внутреннего трения грунта соответственно после их уменьшения относительно фактических 


значений; К, — коэффициент снижения сдвиговой прочности. Величина А соответствующая предельному 


5АМ 2? 
состоянию системы, определяет нижнюю границу прочностных параметров материала. 

2. ГЕМ (Пий едаШбтат тефо4) — метод предельной кинетики, базирующийся на принципе 
Даламбера [8, 9]. Ориентирован на анализ динамической устойчивости сочлененных массивов горных пород. 

3. ТЬЕМ (фи 1Тауег @етепе тефо4 ог ЕЕМ) — метод тонкослойных конечных элементов [10], в котором 
упругопластические элементы тонкого слоя используются для моделирования поведения кинематически 
нестабильных структур. 

Анализ результатов, полученных с помощью методов 5КМ, ГЕМ и ТГЕМ, показал, что в настоящее время 
отсутствует единая концепция математического моделирования поведения структурно нестабильных 
геотехнических систем при нестационарном внешнем воздействии. Это обусловливает актуальность разработки 
методики динамического анализа систем типа «основание — ослабленный слой — блок» в конечно-элементной 
постановке с использованием нового подхода к моделированию плоскостей скольжения. 

Материалы и методы. Уравнение движения механической системы в конечно-элементной формулировке 
представим в виде [11]: 


[МО ЦИН КИ [К ЦИ} =ЗЕз+4Е@)}, (1) 


где [М], [С], [К] — матрицы масс, демпфирования и жесткости ансамбля конечных элементов соответственно; 
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{/"(Е, а, {И} — векторы-столбцы соответственно узловых ускорений, скоростей, перемещений; {АР}, 
Е (Г)} — векторы-столбцы заданных статических и динамических нагрузок соответственно в момент времени #. 


В дальнейшем полагаем, что матрицы [М] и [К] согласованные. 


Для численного интегрирования уравнения (1) используем метод Ньюмарка [12]. При этом шаг 
интегрирования по временной оси ДЕ назначаем так, чтобы с достаточной точностью учитывались вклады 


физически значимых собственных пар. В дальнейшем рассмотрим кинематические способы возбуждения 
колебаний, заданные с помощью либо модельной сейсмограммы \{И’(Р}, либо с помощью модельной 
акселелограммы {И”"(Е)}. При таком способе задания динамического воздействия второе слагаемое правой части 
уравнения (1) будет равно нулю: {А (#)1=0. 

Рассмотрим способ возбуждения механических колебаний посредством модельной сейсмограммы. 
Функцию И’ (ГР) представим в виде [13]: 

Й’ (Г) = А!е`*' зт (04), (2) 

где А — начальная амплитуда; у — коэффициент затухания; 9 — угловая частота внешнего воздействия. На 


рис. 1. приведен график функции И” (2) для значений: А = 0,01553 м; у =0,7143; @=5с (. 


И’ (0,м | 
0,006 





Рис. 1. График модельной сейсмограммы 


Результаты исследования. В качестве первого модельного примера рассмотрим задачу о вынужденных 
колебаниях склона, расположенного на основании (рис. 2). Граничные условия задачи показаны на рис. 3, где 
буквой 5 обозначена исследуемая точка. 


уз Ослабленный 





Рис. 2. Геометрия склона Рис. 3. Расчетная схема склона 


Соотношения геометрических параметров склона и основания (рис. 2) представлены в таблице 1. 
Таблица 1 
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Механические характеристики материала склона и основания: модуль деформации Е = 21 МПа; 
— : С. = -— = : 
коэффициент Пуассона у = 0,3; удельный вес у = 1702 кг/м”; сцепление с = 45 кПа; угол внутреннего трения 


ф = 15°. 
Для моделирования склона и основания используем объемные конечные элементы ЗОГ.ТО45 программного 


комплекса АМЗУ$ Месвашса|. Конечно-элементная модель для варианта с параметрами: /, = 2Н,, Н, = 10 м, 


отнесенная к глобальной декартовой системе координат, приведена на рис. 4. 

Конечно-элементная сетка построена таким образом, что на поверхности контакта смежные узлы 
основания и склона имеют одинаковые координаты, но разные номера. Это принято для того, чтобы разместить в 
этом месте ослабленный слой. Кинематическое воздействие в форме модельной сейсмограммы (2) задаем на 


каждом шаге интегрирования {, в виде узловых перемещений И’(!,) на торцовых поверхностях модели с 
параметрами: Х=Ои Х=А-+/; +[.. 
Моделирование ослабленного слоя (рис. 4) выполняем с помощью упруго-вязких комбинированных 


конечных элементов СОМВГМ14 [14]. Двухузловой элемент СОМВГМ14, состоящий из пружины жесткостью К и 


демпфера жидкостного трения с коэффициентом демпфирования с,, приведен на рис. 5. В рассматриваемом 


случае этот элемент работает только на растяжение-сжатие. 





Х 


Рис. 4. Конечно-элементная модель Рис. 5. Комбинированный конечный 


системы «основание — склон» элемент СОМВПУ14 


В каждом узле контактной поверхности (рис. 4) вдоль глобальных осей Х, У, ГД вводим элементы 
СОМВГУ14. Параметры комбинированных элементов: 
к, =30 кН/м; А, =, =9,44-10 ’ кН/м; с, = 0,5. 


В этом примере и далее вводим допущение о естественном недеформированном состоянии системы 
«основание — ослабленный слой — склон». Для вычислений используем нелинейный решатель комплекса АМЗУЗ 
Месваптса|. 

Результаты конечно-элементного моделирования в виде визуализации деформированного состояния 
системы «основание — склон» с учетом максимального горизонтального смещения и распределения амплитудных 
горизонтальных перемещений И’. (Г) показаны на рис. 6 и 7. Шаг интегрирования уравнения (1) ДЕ = 0,01 с. Как 


видно, введение ЗО) упруго-вязких элементов позволяет смоделировать эффект кинематической нестабильности 
механической системы «основание — ослабленный слой — склон» при кинематическом способе возбуждения 
колебаний. 

й’. (1), м 


-002535 
.004178 
.005822 
. 007465 
.009109 
.010752 
-012395 
-014039 
.0156582 
- 017326 





Рис. 6. Визуализация 


и Рис. 7. Распределение перемещений И’. (1) 


относительно основания 
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Амплитудное значение перемещения в точке 5 составило й = 1,/ см. Для варианта склона 


шах 


[, =Н, (рис.2) И. = 1,1 см. Графики колебаний основания и склона в исследуемой точке 5 (рис. 3) в 


тах 


направлении оси Х приведены на рис. 8. 
-- = 
И, И (0), м 


























0,02 
0,01 
0 
-0,01 
9 Е 
с 


Рис. 8. Графики колебаний в точке 5 основания И”. (Г) исклона 7, (1) 


при кинематическом возбуждении с помощью модельной сейсмограммы 


На основании приведенных графиков видно, что начиная с момента времени Е > 1,5 с наблюдается 
рассогласование колебаний основания и склона. 
Рассмотрим поведение системы «основание — ослабленный слой — склон» (рис. 3) при возбуждении 


колебаний с помощью модельной акселерограммы 7" (1). С этой целью дважды продифференцируем 
выражение (2). В результате получим: 
Й”" (Е) = Ае*' [ у’ зш(0) 1-2 7с03(01)01—зт(01)0“1- (3) 
—2 узш(01)+2с0$(01)6]. 
График функции (3) для параметров: А = 0,01553 м; у =0,7143; 90 =5с ' показан на рис. 9. 


-0,2 
Рис. 9. График модельной акселелограммы 
Кинематическое воздействие в форме модельной акселерограммы (3) по аналогии с сейсмограммой (2) 
задаем на каждом шаге интегрирования #1, в виде узловых ускорений Й”"(#.) на торцовых поверхностях модели с 
параметрами: Х=0 и Х = +1; +1, . На рис. 10 представлены графики колебаний в исследуемой точке 5 (рис. 3) 


при кинематическом воздействии в форме модельной акселерограммы. Сравнивая графики колебаний, 
приведенные на рис. 8 и 10, устанавливаем, что они практически совпадают. Это свидетельствует о корректности 
разработанной конечно-элементной модели, позволяющей описывать поведение системы «основание — 
ослабленный слой — склон» при различных способах нестационарного кинематического воздействия. 
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Рис. 10. Графики колебаний в точке 5 основания И”.. (Г) исклона 7, (1) 


при кинематическом возбуждении колебаний с помощью модельной акселелограммы 


В качестве второго модельного примера рассмотрим задачу о вынужденных колебаниях склона с 
кинематически нестабильным клиновидным включением (рис. 11). В связи с симметрией конфигурации, в 
расчетной схеме учитываем только 1/2 часть склона и включения. Граничные условия для принятой расчетной 
схемы показаны на рис. 12. Здесь буквой 5 обозначена исследуемая точка, принадлежащая одновременно 


основанию склона и клиновидному включению. 
у 







Ослабленный 
слой 
Направление 
эх скольжения 
Й’ =0 
р й=0 52 
У 
Рис. 11. Схема склона с клиновидным Рис. 12. Расчетная схема для задачи «склон — 
включением клиновидное включение» 


Конечно-элементная модель склона и клиновидного включения показана на рис. 13. Как и в предыдущем 
примере, в данном случае используем элементы 5ОТ.1245 и СОМВГ\М14 с теми же характеристиками материала. 


ь 


а) 6) в) 


Рис. 13. Конечно-элементная модель: а — склон; 6 — клиновидное включение; 





в — склон с клиновидным включением 


Функция, описывающая модельную акселелограмму, имеет вид: 


3 

— ки 

Й’ (1) =А-со$ бЕ, С 

ти сз 

где А — амплитуда ускорения; 9 — частота внешнего воздействия. На рис. 14 приведен график И’ (Г) при А = БК 
2,5 м/с’, 09 = 2Гц. Значения ускорения И (&.) на 1-том шаге интегрирования уравнения движения (1) 7 


прикладываем к узлам поверхности модели с координатой Х = 0 (рис. 12). 
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ИР. (В), м/с? 





0 0,05 0,1 0,15 0,2 


Рис. 14. График модельной акселелограммы И (1 ) 


Результат моделирования в виде распределения амплитудных значений перемещений И’. (Г) показан на 
рис. 15. Шаг интегрирования ЛЕ = 0,01 с. Графики колебаний основания склона и клиновидного включения в 


исследуемой точке 5 (рис. 12) в направлении оси Х приведены на рис. 16. 
Как видно из рис. 15 при заданном кинематическом воздействии происходит разрыв сплошности массива 
склона по ослабленному слою и клиновидное включение смещается относительно основания склона вдоль оси Х. 


Й’, (#),м 


—-02252- 
—.0198'74 
—-017224 
—-014573 
—-011923 
—.009273 
—.005622 
—.003972 
—.001321 
„001329 


У 


5 


Рис. 15. Распределения И’. (1) в 1/2 части склона с клиновидным включением 
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Рис. 16. Графики колебаний в точке 5 основания склона И/^ (1) и клиновидного включения И’ (1) при кинематическом 
5 5 


возбуждении колебаний с помощью модельной акселелограммы 


«Дрейф» И”. (г) на рис. 16 обусловлен тем, что данная конечно-элементная модель не имеет связей, 


препятствующих смещениям вдоль оси Х. Как показано в [11], решить проблему «дрейфа» можно путем вычитания из 
значений перемещений И’. (Г) и И, (1) перемещения основания склона, которое представляет собой смещение «как 


жесткое целое». Отметим, что полученные амплитудные значения перемещений И’, (1) и И’, (1) позволяют оценить 
динамические параметры системы «склон — ослабленный слой — клиновидное включение». 


Заключение. Разработана и верифицирована конечно-элементная модель для исследования 
динамического поведения кинематически нестабильных склонов в трехмерной постановке с учетом физической 


Гайджуров П. П. и др. Численное моделирование поведения кинематически нестабильных склонов при динамических воздействиях 





нелинейности материала. 
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